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ABSTRACT 



Aims. We aim to study dust properties of massive star forming regions in the outer Galaxy, in a direction opposite to the Galactic 
center 

Methods. We present observations of six outer Galaxy point sources IRAS 01045+6505, 01420+6401, 05271+3059, 05345+3556, 
20222+3541 and 20406+4555, taken with the Submillimeter Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell 
Telescope (JCMT) at 450 and 850 /jm. Single temperature greybody models are fitted to the Spectral Energy Distribution of the 
detected sub-mm cores to derive dust temperature, dust emissivity index and optical depth at 250 /im. The observed radial intensity 
profiles of the sub-mm cores were fitted with power laws to derive the indices describing the density distribution. 
Results. At a resolution of 15" all six IRAS point sources show multiple emission peaks. Only four out of fourteen detected sub- 
mm cores show associated mid-infrared emission. For the sub-mm cores we derive dust temperatures of 32±5 K and dust emissivity 
indices between 0.9 and 2.5. The density profiles of the sub-mm cores can be fitted by a single power law distribution with indices 
-1.5±0.3, with most cores showing an index of -1.5. This is consistent with most observations of massive star forming regions and 
supports predictions of models of star formation which consider non-thermal support against gravitational collapse. 
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1. Introduction 

With the development of high angular resolution millimeter 
(mm) and sub-millimeter (sub-mm) instrumentation, observa- 
tions studying the influence of different physical environment 
on the formation of molecular clouds subsequently developing 
into protostellar cores has been receiving increased attention. 
While the inner Galaxy has been extensively studied, lack of 
instruments with high angular resolution has limited the study 
of sources in the outer Galactic disk which are in a poten- 
tially different physical and chemical environment. Existing ob- 
servations suggest that in the outer Galacti c disk the molec- 
ular clouds are more sparsely distributed ( IWouterloot et al.L 
Il990l). the diffuse galactic inters tellar radiation field is weaker 
(Cox & Mezger, 1989; Bloemen, 1985), the m etalUcity is lower 



(Shaver et al., 1983; Fich & Silkey, 1991; Rudolph et a l.Ll2006[) 
and t he cosmic-ray flux density is smaller ( Bloeme n et al.L 
Il984 l). compared to the solar neighbourhood. Although the effi- 
ciency at which molecular clouds form from atomic gas is much 
lower in the far outer Galaxy, the mass spectrum of molecu- 
lar clumps i s not significantly different from the inner Galaxy 
dSnell et al.L |2002|) . The spectra however show a tendency to- 
wards steepening in the far outer Galaxy. The star formation 
activity in the outer Galaxy clouds is also similar to the in- 
ner Galaxy clouds. Thus the observed low global star forma- 
tion rate in the outer Galaxy is a consequence of the inefficiency 
of the formation of molecular clouds in this region and not an 
inefficiency in the star formation process within these clouds 
(Snell et al., 2002). The outer Galaxy star forming regions could 



Send offprint requests to: bhaswati@tifrres.in 



thus be better examples of stars and/or stellar clusters forming 
in isolation in contrast to the crowded fields within the inner 
Galaxy. It is thus important to study the point sources identified 
by IRAS in the outer Galactic disk at high angular resolutions 
to probe their multiplicity, dust properties and kinematics. As a 
first step towards these goals we present here sub-mm continuum 
mapping of selected outer Galaxy IRAS sources. 

We have selected a sample of outer Galaxy (with 
Galactocentric distances greater than 10 kpc) sources from the 
IRAS point source catalog (PSC) based on the following cri- 
teria: Fioojum > 100 Jy, F,oo^m >F60jum, 30°< / <330°, 
-5°< b <5°. The Galactic longitude criterion ensures that 
source confusion along the line of sight faced towards the in- 
ner Galaxy, is avoided. The sample consists of 17 sources, all 
of which are bona-fide point sources according to the IRAS 
PSC. For majority of these sources no continuum maps other 
than the IRAS images are available. Here we present sub- 
mm (450 and 850 yum) observations of six sources in our 
sample, viz., IRAS 01045+6505, IRAS 01420+6401, IRAS 
05271+3059, 05345+3556, 20222+3541, and 20406+4555. For 
the sources IRAS 01045+6505 and IRAS 01420+6401 we have 
only 850 jum observations, hence in this paper, we present only 
the basic results for these two sources and discuss the other four 
sources in detail. 

The purpose of these observations was primarily to obtain 
a higher angular resolution view in order to detect the sub- 
structures within these sources. Using these observations we also 
derive the dust properties of these regions in order to improve our 
understanding of the true nature of sources located in the outer 
Galaxy environments. 
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2. Observations and data reduction 

We present here results of observations made between 
September 2004 and January 2005, using the Submillimeter 
Common-User Bolometer Array (SCUBA) on the James Clerk 
Maxwell Telescope jHoUand et all L1999). The SCUBA aiTay 
covers a hexagonal 2'.5 field of view with 9 1 and 37 pixels at 450 
and 850 fim respectively. Maps are fully sampled using the "jig- 
gle" mode, in which the telescope is moved around a 64-position 
pattern by the secondary mirror in order to fully sample the sky 
with the wider intrinsic beam spacing of the array. A chop throw 
of 120" in azimuth was used for all observations. 

The data were r educed using the SCU B A Use r Reduction 
Facility (SURF; iJenness & Economoui [1999*) and the 
STARLINK imaging software follo wing the methods described 
in the SCUBA mapping cookbook dSandell. Jessop & Jenness L 
l200lh . The steps of reduction included flat-fielding, extinc- 
tion correction, sky-subtraction and calibration of images in 
Jy beam"'. The maps were calibrated using observations of 
primary calibrators like Uranus, CRL 618 and CRL 2688 
observed close to the time of our observations. 

All the maps were converted to FITS -files and exported to 
MIRIAD (Sault et al., 1995) for further analysis. In order to cor- 
rect for the error lobe contribution, we deconvolved all the maps 
using CLEAN and a symmetry-sized model beam. The model 
beam, derived from the 450 and 850 yum observations of Uranus 
obtained during the same night and under similar conditions, and 
consistent with repeated observations during the whole observ- 
ing run, is composed of three symmetric Gaussians. At 850 fim 
it has HPBWs of 15'.'2, 58'.'6 and 140", with ampHtudes 0.985, 
0.013 and 0.002, respectively. At 450 jum it has HPBWs of 8'.'5, 
33" and 140" with amplitudes 0.942, 0.055 and 0.003 respec- 
tively. In order to obtain images with good S/N we have restored 
maps at both 450 and 850 yum to a common resolution of 15". 

3. Results 

Figure [1] shows the 450 and 850 jim dust continuum maps of 
the six outer Galaxy IRAS point sources. We find that all the 
six sources are resolved into at least two to three emission peaks 
in the sub-mm images. Comparison with the MSX point source 
catalog suggests that most of the sub-mm cores do not have 
any associated MSX source. Given the 4" positional accuracy 
of the SCUBA maps (see below), we consider a sub-mm core 
to be associated with an MSX source only if the MSX source 
is within 5" of the position of the sub-mm core. Given the 
large distances to the sources, displacement of the young star 
from its core due to proper motion outside of this association 
range is rather unlikely. Our interpretation is that the lack of 
MSX counterparts for the sub-mm sources is due to differences 
in the evolutionary stages. Based on this criterion we find that 
only 05271-SMMl, 05271-SMM2 and 01045-SMMl have as- 
sociated MSX sources. Only IRAS 01045-1-6505 was observed 
using IRAC and MIPS instruments on Spitzer. Preliminary re- 
sults from Spitzer show that bofli 01045-SMMl and 01045- 
SMM2 have sources at 24 yum associated with them, with the IR 
source at the position of 01045-SMMl being partially saturated. 
This indicates that 01045-SMMl is a reasonably hot source 
with a dust temperature of 40 K or more. IRAS 20406+4555, 
01420-H6401 and 01045+6505 have associated radio continuum 
sources. For IRAS 01420+6401 and 01045+6505 the radio con- 
tinuum positions are b ased on high resolution VLA observations 
( lRudolphetaIlll996h . wher eas for IRAS 20406+ 4555 the radio 
data is from the GB6 catalog jGregorv et al.L[T996l) . The nominal 



IRAS position for 05345+3556 is offset with respect to the ridge 
detected in the sub-mm. We have checked the ISSA images for 
this source and find that at the position given by the IRAS PSC 
there is no prominent peak. In Sect. |4] we discuss the sources 
individually in the light of the available limited complementary 
observations. 

Given the extended nature of these sources and the 2' chop 
throw, the observed off positions will not have been completely 
free of emission, and this will have been subtracted from the 
fluxes and as a background from the images. However, as we are 
more interested in the warmer dust corresponding to the warm 
cores seen by IRAS, rather than the cold diffuse component, 
the background-subtracted flux measurements are appropriate. 
In order to derive the parameters of the sub-mm cores detected 
in the IRAS sources, we have fitted an elliptical Gaussian for 
the sub-mm source, and a background plane for the surrounding 
cloud using the task IMFIT in MIRIAD. The fit to the broader 
component is mainly to provide a good subtraction of the ex- 
tended emission, and is not to estimate the flux density of the 
surrounding cloud. Table [Upresents results of the fits to the de- 
tected sub-mm sources. We estimate the uncertainty in determin- 
ing the positions of the sub-mm sources based on the Gaussian 
fitting to be 1". Combining this with the pointing accuracy of 
the observations (< 3") we estimate the final astrometric uncer- 
tainty of the positions of the sub-mm cores to be ~ 4". The sizes 
are derived from the Gaussian fits to the 450 fim data and the 
flux densities at 850 yum are derived by constraining the sizes of 
the Gaussian fits to be the same as the sizes at 450 yum. For the 
sub-mm cores not clearly detected at 450 yum, we have derived 
the sizes from the 850 yum map and integrated the 450 fim fluxes 
over an aperture of same size. 

4. Notes on the sources 

Here we present details of existing observations of the four 
sources IRAS 05271+3059, 05345+3556, 20222+3541 and 
20406+4555. 

IRAS 05271+3059: This source, WB89 651, is part of the 
outer Galaxy catalog compiled by Wouterloot & Brand ( 198^. 
Based on the velocity information from the molecular line obser- 
vations the distance to the source is estimated to be 16.5 kpc. The 
luminosity of the source estima ted from the 60 and 10 fim IRAS 
flux es is 4.2 10^ Lp. C S(2-1) (Bronfman etaD, [11%) and H2O 
(ValdettaroetaUl2001h maser emission was found to be associ- 
ated with this sour ce, although no methanol maser was detected 
(ISlvshet all [19991) . There is no radio continuum source known 
to be associated with IRAS 05271+3058. We find that both sub- 
mm sources are associated with MSX point sources (Fig.[T]i- 

IRAS 05345+3556: Also known as WB89 673, this source 
is at a distan ce of 14 kpc as derived from molecular line ob- 
servations by IWouterloot & Brandl (Il989b and has a luminos- 
ity of 1.7 10^ Lq. The radio continuum sour ce detected in 
NRAO VLA Sky Survey (NVSS) at 1.4 GHz dCondon et all 
1998) is coincident with the nominal position of the infrared 
peak according to the IRAS Point Source Catalog. As noted in 
Sect. [3] the dust emission detected in the SCUBA maps is along 
a ridge which is offset from this nominal infrared peak posi- 
tion. Simila r to IRAS 05271 +3059, CS (2-1) ( Bronfman et all 
1996), SiO jHariu et al.L [19981) and H2O jValdettaro et al..,.200lh 
maser emissions are found t o be associated wi th this source and 
no methanol was detected ( ISlvsh et all Il999l) . No MSX point 
sources were found within a search radius of 20" around the 
sources SMM-1 and SMM-2, only SMM-3 has a nearby MSX 
source. 
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Fig. 1. Contours of 450 & 850 jjm emission for the outer Galaxy sources. '+' and * indicate positions of MSX and IRAS point 
sources respectively in the mapped regions. Filled triangles show locations of radio continuum sources. The peak flux densities 
(rms) are as foUows: (1) 05271+3059: Fgso = 0.6 (0.02) Jy/beam, F4go = 4.7 (0.03) Jy/beam; (2) 05345+3556: Fgjo = 1-0 (0.03) 
Jy/beam, F450 = 8.0 (0.04) Jy/beam; (3) 20222+3541: Fsso = 0.5 (0.02) Jy/beam, F45() = 3.2 (0.02) Jy/beam; (4) 20406+4555: Fsso 
= 0.5 (0.02) Jy/beam F450 = 1.8 (0.01) Jy/beam; (5) 01420+6401: Fsso = 0.4 (0.02); (6) 01045+6505: Fgso = 1-2 (0.03) Jy/beam. 11 
equally spaced logarithmic contours are drawn, with the lowest level at 2cr and the highest level at 95% of the peak value. 



IRAS 20222+3541: iBeck & BroschI (Il988l) suggested fliat 
IRAS 20222+3541, also known as MCG 06 -45-001, is a lo- 
cal group galaxy. In a subsequent paper. Shor e & Sagd d 19901) 
used CO observations to conclusively show that rather than be- 
ing an external galaxy, this IRAS source is a Galactic star form- 
ing region wi th an associate d H n r egion. Based on the CO ob- 
servations of Shore & Sage (1990) we estimate the kinematic 
distance to the source to be 2.2 kpc for a visr of ~7.4 km/s. 
This is the nearest of the sources in our sample. The luminos- 
ity of the source is 1.5 10^ L©. Only one point source from the 
MSX catalog is found in the mapped region although the im- 
ages of almost all MSX bands show an extended emission co- 
incident with the IRAS source. The only identified MSX point 
source in this region i s offse t from SMM-1 by ~ 25". In con- 
trast to IShore & Sagd d 19900 we find that the only radio con- 
tinuum sour ce detected in t he NR AO VLA Sky Survey (NVSS) 
at 1.4 GHz ( Condon et"aLlll998l) is off'set by ~ 2' with respect 
to the dust continuum peaks seen in IRAS, MSX and SCUBA 
datasets. In 



(Il989h . is located at a distance of 7.8 kpc. We have confirmed 
from our recently concluded and yet unpublished CO observa- 
tions, as well as from the '^CO J=l-0 observations of IWu et akl 
(2001) that only the velocity component at -49 kms"' is associ- 
ated with this source and this corresponds to the aforementioned 
distance. The luminosity of the source based on the IRAS fluxes 
is 2.6 10"^ Lq. There are no MSX point sources in the MSX6C 
associated with this region, although both in bands A and E the 
source is clearly detected. An H ii reg ion from the GB6 catalog 
of radio sources (^ Gregory et all 1 19961) is centerd at the position 
of SMM-1, while the IRAS source is located closer to SMM-2. 
The comparatively compact appearance of SMM-2 & 3 together 
with the H 11 region might suggest that SMM-2/3 contain proto- 
stars, formation of which is triggered by the H 11 region. However 
so far n o NH3 and H2O mase r emission were detected from this 
source dMolinari et al.Lll996h . 



IRAS 20406+4555: This source, als o known as WB89 4 
from the catalog of outer Galaxy sources bv lWouterloot & Brandl 
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Table 1. Observed properties of the sub-mm sources identified in the SCUBA maps 



Submillimeter 
Source 


Q'(2000) 


5(2000) 


D 

kpc 


e„ X Bi," 

("X") 


S450 

Jy 


Sgso 

Jy 


IR 

Assoc. 


Radio 
Assoc. 


05271+3059 SMM-1 
05271+3059 SMM-2 


05:30:20.97 
05:30:20.28 


31:01:25.8 
31:00:50.3 


16.5 
16.5 


12.9 X 12.8 
15.1 X 12.7 


7.92 
2.26 


1.08 
0.24 


MSX 
MSX 




05345+3556 SMM-1 
05345+3556 SMM-2 
05345+3556 SMM-3 


05:38:00.41 
05:38:01.29 
05:38:00.99 


35:59:00.2 
35:58:21.2 
35:57:45.7 


14.0 
14.0 
14.0 


19.2 X 14.1 
13.9 X 9.2 
17.7 X 14.5 


14.18 
12.74 

b 


1.65 
1.14 
0.80 






20222+3541 SMM-1 
20222+3541 SMM-2 


20:24:10.69 
20:24:11.85 


35:51:10.2 
35:50:34.8 


2.2 
2.2 


20.4 X 20.0 
15.0 X 13.2 


6.06 
5.95 


0.33 
0.51 






20406+4555 SMM-1 
20406+4555 SMM-2 
20406+4555 SMM-3 


20:42:21.75 
20:42:21.30 
20:42:19.33 


46:06:21.1 
46:05:42.8 
46:05:38.3 


7.8 
7.8 
7.8 


15.4 X 11.0 
12.7 X 12.6 
15.7 X 11.3 


1.79 
3.02 
2.17 


0.17 
0.30 
0.24 




GB6 


01420+6401-SMMl 
01420+6401-SMM2 


01:45:39.09 
01:45:40.76 


64:16:15.9 
64:15:57.3 


10.1 
10.1 


14.8x10.8 
16.5x16.3 




0.44 
0.71 




GB6,VLA 


01045+6505-SMMl 
01045+6505-SMM2 


01:07:50.99 
01:07:54.99 


65:21:23.2 
65:20:31.7 


11.3 
11.3 


12.7 X 11.6 
16.7 X 16.5 




1.81 
0.46 


MSX,Spitzer 
Spitzer 


GB6,VLA 



" Deconvolved Major (9^) and Minor (0^) axes derived from elliptical Gaussian fits. 

' Not detected at 450 fim. S450 obtained by integrating at the position of the 850 jjm peak is 5.29 Jy. 



5. Spectral Energy Distributions of the sub-mm 
cores 



Several authors ('Sridharan et alj 120021; iBeuther et all |2002|) 
have shown that two temperature components are present in the 
SED of luminous YSOs: a compact hot component which dom- 
inates the IRAS 12 yum and 25 yum fluxes, and a more extended 
component arising from colder gas which dominates the IRAS 
60 yum and 100 fim fluxes. Here we analyze the spectral en- 
ergy distributions (SED) of the sub-mm cores detected at both 
SCUBA wavelengths, between 60 and 850 i-im in order to derive 
order of magnitude estimates of the dust temperature and mass 
of the cold component. Since the sub-mm cores are not resolved 
in the mid and far-infrared, we apportion the 60 and 100 yum flux 
densities for the IRAS sources to the sub-mm cores in the same 
ratio as the observed 450 jim fluxes. 

The greybody functions have the form (e.g., 

iDent. M atthews, & Ward-T hompsoa . 1 998,1 . 



(1) 



where Fy is the flux measured at a frequency v, Q is the solid 
angle subtended by the cloud, By(Tdust) is the Planck function 
evaluated at the dust temperature Tdand v, and Ty is the optical 
depth at v. The optical depth Ty scales as tqv^, where /3 is the 
grain emissivity index and tq is the optical depth at a reference 
frequency (vq). 

Fitting was performed such that was minimized, with 
defined as. 



X 



V, model 



v,ohs 



(2) 



This definition o f gives equal we ighting to the different 
wavelength regimes dHunter et aI.L 12000'). We have constrained 
the angular sizes of the sub-mm cores to those measured in 
the SCUBA observations (Table |2]l. Three free parameters were 
used: the temperature, the optical depth at 250 yum and the grain 
emissivity index (/?). 

Figure |2] shows the greybody models fitted to the SEDs 
and Table |2] also presents the parameters like mass, luminos- 
ity etc. calculated for the best fit model. The mass has been es- 
timated from the fitted values of T250, Ta and /3, adopting the 
"Hildebrand" mass opacit y, Kq, defined at 250 yum (1200 GHz), 
i.e. /fi200GHz =0.1 cm^g"' (iHildebrandi 1 1 98 3l) and a gas-to-dust 
ratio of 100. The main uncertainty in calculating the dust masses 
here lie in the use of inner Galaxy dust properties as well as 
gas-to-dust ratio for these outer Galaxy sources. Observationally 
there is evidence of re duced metallicity towards the outer Galaxy 
(iRudoIph et aI.L 12006") and the dust-to-gas ratio depends almost 
linearly on metallicity (Boselli et al., 2002). The assumed gas- 
to-dust ratio of 100 does not take the metallicity gradient into 
account. In addition, the calculated dust mass opacities depend 
on the grain composition and even for inner Galaxy sources the 
commonly used values differ by as much as factors of 3. Overall, 
we estimate that the calculated mass can deviate from the real 
value by up to a factor of 10. For the sub-mm sources in IRAS 
01045+6505 and 01420+6401 in the absence of the 450 //m 
fluxes greybody SED fitting was not possible. We have estimated 
the mass of these sub-mm cores assuming to be 32 K, average 
value of Td found from the SED fits and (3 - 2. As discussed in 
Sect. |6]a choice of p - 2 can be reasonably justified based on 
both theoretical and observational results. 
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Table 2. Parameters of the SED Models for the sub-mm cores 



Source 


Dia 


T 


■r250 


P 


L 


M 




pc 


K 






Lo 


Mo 


05271-SMMl 


1.03 


36.9 


0.01 


1.2 


42272 


711 


05271-SMM2 


1.10 


28.6 


0.01 


1.9 


15184 


814 


05345-SMMl 


1.12 


30.7 


0.01 


0.9 


13192 


595 


05345-SMM2 


0.11 


28.6 


0.03 


1.2 


12270 


778 


05345-SMM3 


1.09 


28.6 


0.01 


1.0 


6477 


398 


20222-SMMl 


0.22 


26.6 


0.02 


2.5 


1091 


86 


20222-SMM2 


0.15 


28.6 


0.03 


2.2 


1083 


61 


2U406-SMM1 


0.49 


28.6 


0.01 


2.4 


6448 


320 


20406-SMM2 


0.48 


30.7 


0.02 


2.4 


12463 


422 


20406-SMM3 


0.50 


32.8 


0.01 


2.4 


15387 


334 


01420-SMMl 




32° 








1041* 


01420-SMM2 




32" 








350* 


01045-SMMl 




40" 








267* 


01045-SMM2 




32" 








432* 



" Assumed T^. 

* Mass calculated from Sgso assuming the given in the table and 
dust emissivity index = 2. 



Except for the sub-mm cores in IRAS 05345+3556, the grey- 
body models produce reasonable fits to the SEDs of the sub-mm 
cores in IRAS 05271+3059, 20222-H3541 and 20406+4555. The 
SED models for the cores in IRAS 05345+3556 do not repro- 
duce the 450 im\ fluxes. A somewhat less constrained two dust 
temperature greybody model fit to the SED of the sub-mm cores 
in IRAS 05345+3556 yields better fits. This indicates that ad- 
ditional observations between 100 and 450 jum are essential to 
fully constrain the fits to the SEDs. 



6. Dust Properties from the SED fits 

The temperatures of the cold components for all the sub-mm 
cores are ~ 32+5 K and the dust emissivity indices vary between 
0.9 and 2.5. While the sub-mm cores within IRAS 05271+3059 
show very different dust emissivity indices (/?; Table|2]i, the three 
sub-mm cores in IRAS 20406+4555 show almost identical val- 
ues of p. 

The dust grain emissivity index is crucial in determining 
many source parameters, including the dust temperature, as 
well as providing information about the gr ain structure and 
the dielectric characteristics of the grains (ISchwartzL [1981 
iHildebrandl Il983l) . Knowledge of y6 improves estimates of the 
dust opacity and subsequently of the dust mass of the star 
forming clouds. Various estimates of fi from both laboratory 
and theoretical models have been prop osed and so far there 
is little consensus between these studies. iDraine & Led d 19841) 
derived /3 ~ 2 for 40< A < 1000 //m, using the avail- 
able laboratory data to measure optical constants of a mix- 
ture of naked graphite and silicate grains. A dust grain emis- 
sivity index of ~ 2 has als o been supported by the calcula- 
tions of [K nacke & Thomso nI (ll973l):lM athis & Whiffen ( 1989); 
iKnipel & Sieb enmorgeni (|1994|) . iMivake & Nakagawai 619931) 
howe ver d erived values of /3 in excess of 2 for large grain sizes, 
while Wright ( 1987) calculated 0.6< j8 < 1.4 for fractal grains. 
lAannestad (.1975.) predicted f3 up to 3.5 for olivine, fused quartz. 



and lu nar rock grains covered with ice mantles. iMennella et al.l 
(Il995h concluded that more evolved cores have a higher emis- 
sivity index (J3). 

Observationally there is an equally large scatter in the de- 
rived values of dust emissivity index, /3, as in theoretical stud- 
ies. There are several instances of obs ervations of circum- 
stellar disks around young PMS stars (Beckwith & Sargenj, 
1991; Sandell & Weintraub, 2001; Natta et al., 2004), low mass 
protostars dHogerheijde & Sandell L 12 000; Stark et al., 2004) as 
well a s high mass protostars jSandelll i2000; .Sandell & Sieversl 
120041) for which /3 < 2 were derived. L arger values of 13 , 
i.e., in excess of 2 hav e also been reported (Kuan et al.l 119961 ; 
Goldsmith, Be rgin. & L is. 1997; Lis & Menten, 1998). We note 
here that the derived values of /3 is always to some extent af- 
fected by the temperature and density variations as long as there 
is some optically thick portion of dust and detailed radiation 
transfer mode ls can achieve a good fit to observations using 
yS = 2 (Mueller et al.L 120021; IS hirlev et all 12002). In addition, 
uncertainties and artifacts due to ignored error-lobe pickup, cal- 
ibration uncertainties, optical thickness of H ii regions, unequal 
beamsizes at the two wavelengths under consideration etc. have 
often contributed to observationally determined higher values of 
P- 

The values of /3 derived from the SCUBA observations of the 
outer Galaxy sources lie within the limits of the different values 
seen in the inner Galaxy sources. Based on the present estimates 
we conclude that the overall dust properties shown by the outer 
Galaxy sources, do not differ significantly from the properties 
of the inner Galaxy star forming regions, and show no indica- 
tion of the substantially higher values of /3 as found in a few 
extreme regions in the Galactic Center environment. The proto- 
stellar sources observed here are rather distant, hence the dust 
emission is dominated by the extended dense cloud envelopes 
surrounding the young high mass protostellar objects, rather than 
being dominated by the emission from the accretion disk, which 
are always found to show a lower apparent /3. 
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Fig. 2. Single temperature greybody fits to the Spectral Energy 
Distribution (SED) of sub-mm cores. Fitted parameters are pre- 
sented in Table 12] 



7. Radial Profiles 

7.1. Fits to the radial intensity profiles 

Figure |3] shows the radial intensity profiles at 850 fim of all the 
sub-mm cores with the exception of 05345-SMM3 and 20406- 
SMM3. These radial profiles have been derived by averaging all 
datapoints within a 2" wide circular annulus at a radius r from 
the center of the core, with r ranging from to 30". The obser- 
vations were done with a chop throw of 2' in the azimuth. Based 
on the derived angular sizes of the sub-mm cores we conclude 
that the effect of chopping on the derived radial intensity profiles 
is negligible. 

We have created model intensity distribution by assuming an 
analytical two-dimensional power-law intensity distribution of 
the form I{r) oc 0^ and convolving it with the synthetic beam pro- 




Radial distance from sub— mm peak ["] 

Fig. 3. Normalized, azimuthally averaged radial profiles for the 
intensities of the sub-mm cores at 850 //m (points). Solid line 
shows the best fit power-law model corresponding to the index 
y. The dashed line shows the model beam profile. The errorbars 
indicate the rms values within the annuli. 



file at 850 jum which we used to deconvolve the intensity maps. 
We derived the model radial intensity profile by azimuthally av- 
eraging exactly in the same manner in which we have created 
the observed intensity profile. We have performed a.x^ fitting of 
a family of such model profiles corresponding to different power 
law indices to the observed model profile. Keeping in mind the 
angular resolution of the 850 images we have restricted the 
fitting procedure to radii 8"< r < 30", so as to have a reliable es- 
timate of the nature of the intensity distribution at spatial scales 
unaffected by the resolution of the observations. Figure[3]shows 
the resultant fits and Table [3]presents the fitted values of j. 

For most of the cores it is possible to derive acceptable fits to 
the radial intensity profiles of the outer envelopes (r > 8") using 
a single power law. The fits, though not constrained by the inner 
regions, are reasonably good at radii <8" as well and we derive 
y - -0.9 + 0.3 for the 8 sub-mm cores. 

7.2. Density Profiles of the sub-mm cores 

Here we consider a model with a spherically symmetric enve- 
lope and use the intensity distribution of the dust emission in the 
SCUBA images at 850 /im to derive a first estimate of the radial 
density profile in the envelope. The density profile is a power 
law, n oc rP, the value of p derived using this method depending 
upon the assumed temperature distribution. 

For optically thin dust emission, assuming the Rayleigh- 
Jeans approximation to be valid, spherical cores and tempera- 
ture and density distribution following a power law, the intensity 
index y depends to first order on the d ensity index p and the 
temperature i ndex q (T cc r*) through jMotte & Andrei 120011 : 
lAdamsL[T99l : 



y^l-p-q 



(3) 



Thus to derive the density index p it is necessary to have 
some knowledge about the temperature distribution. Based on 
the grey-body fitting of the SEDs it is apparent that the sub-mm 
cores are centrally heated (accretion and/or stellar burning) so 
that the temperature decreases with increasing distance from the 
center. Theoretical calculations of massive star forming regions 
suggest that steeper indices for the temperature distribution are 
expected in the inner regions (<2000 AU). However beyond, the 
inner 2000 AU, the temperature distribution flattens aga in and 
runs asymptotically into a r"" "* distribution ( van der Tak et al.L 
2000). For the large distances to our sources, the typical reso- 
lution of the 850 i-im maps are well beyond 20,000 AU. This 



Mookerjea et al.: SCUBA studies of Outer Galaxy IRAS sources 



7 



Table 3. Power law indices for radial intensity & density profile 



Source 


r 


P 


05271-SMMl 


-1.2 


-1.8 


05271-SMM2 


-0.9 


-1.5 


05345-SMMl 


-0.9 


-1.5 


05345-SMM2 


-1.0 


-1.6 


20222-SMMl 


-0.5 


-1.1 


20222-SMM2 


-0.9 


-1.5 


20406-SMMl 


-0.7 


-1.3 


20406-SMM2 


-0.9 


-1.5 



implies that we can safely assume the asymptotic solution of 
q - 0.4 for our sources. Thus, for parts of the sub-mm cores at 
radii larger than 20,000 AU we derive density indices varying 
between -1.1 and -1.8, with most cores showing indices of -1.5 
(Tableig. 

Analytic and numerical studies show that the mass infall rate 
during the protostellar phase is dependent on the radial density 
profile at the onset of collapse, and on the equation of state of 
the material. iMotte & Andrei ( |200ll) studied isolated and clus- 
tered protostellar envelopes, and found that isolated protostel- 
lar envelopes exhibit density profiles like n oc r~^, as predic ted 
by the standard model of star formation bv lShu et~aL I (ll987h .ln 
contrast, protostellar envelopes in clusters are found to be in- 
duced in compact condensations resembling more finite-sized 
Bonnor-Ebert spheres than singular isothermal spheres, which 
suggests that dynamical protostellar models are more appro - 
priate (Mott e & Andri, 1200 ll: IWhitworth & Summersl 1 19851) . 
ISO based mid-infrared absorption studies of low-mass prestel- 
lar cores revealed a flattening of the inner parts of the radial 
profiles as well as a s teepening of the profiles farther outside 
dBacmann et a n i2000l) . In the high mass regime studies of ul- 
tracompact H ii regions and/or hot cores indicate that the overall 
density profiles are more like n oc with -1.0 < p < -1.5 as 
predicted by loga t ropic equations of state dHatchell et alll2000l; 
Ivan der Tak et all 12000 ). Based on the stu dy of a large s ample 
of high mass protosteUar objects (HMPOs) lBeuther et all d2002l) 
derived a much wider range of density indices n ~ ^-i-6±0 5 ^ 
possible physical interpretation of the flattened density profile 
for massive objects could be that in low-mass objects the sup- 
port against gravitational collapse is thermal, while for massive 
objects no nthermal support p lays a major role (Myers & Fuller, 
[1992; McLau ghUn & Pudritzi ll996). 

Based on the luminosity estimates (Table |2| the sources we 
consider here are all intermediate to massive star forming re- 
gions, some of them having Hii regions associated with them 
(Fig.[T]|. The derived value of density indices for the 8 sub-mm 
cores is -1.5 ± 0.3 and is thus consistent with observations of 
massive star forming regions and agree with theoretical models 
with significant non-thermal support against gravitational col- 
lapse. 

8. Summary 

The high angular resolution (15") sub-mm dust continuum maps 
of the outer Galaxy IRAS point sources show that all of them 
contain multiple emission peaks. Only four of the fourteen de- 
tected sub-mm sources have associated mid-infrared emission 
as seen from the MSX point source catalog and/or Spitzer IRAC 
and MIPS observations. SEDs for individual sub-mm cores be- 
tween 60 and 850 yum have been generated by apportioning the 
unresolved IRAS fluxes at 60 and 100 fj.m in the same ratio 



as the observed 450 //m fluxes. Single temperature greybody 
models were fitted to the observed SEDs of the sub-mm cores 
to derive the dust temperature, emissivity index (J3) and opti- 
cal depth at T25o- The derived dust temperature of the sub-mm 
cores is 32+5 K and the emissivity index ranges between 0.9 and 
2.5. The two sub-mm cores in IRAS 05271-1-3059 show largely 
different emissivity indices. In spite of the known radial gradi- 
ents in metallicity within the Milky Way, the dust properties de- 
rived here do not show any significant difference from the inner 
Galaxy star forming regions. 

An analysis of the radial intensity profiles (for r > 8") at 
850 jum for the sub-mm cores shows that it is possible to fit 
the profiles reasonably well using a single power law with in- 
dices -0.9+0.3. Deriving density distributions from the radial 
indices we get density power laws -1.5+0.3, with most cores 
showing an index of -1.5. This density profile is similar to the 
profile observed in many massive star forming regions in the in- 
ner Galaxy and generally agree with models of star formation 
which invoke non-thermal support against gravitation al collapse 
dMvers & FufleA 1 19921; [McLaughlin & Pudritzll 19961) . 

Thus, except for their isolated occurence in the outer Galaxy, 
in most other respects the sources studied here look very simiilar 
to star formation cores in the inner parts of the Galaxy which 
have similar range of masses. 
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